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Abstract. The paper addresses the homogenization of a family of micro- 
models for the flow of a slightly compressible fluid in a poroelastic matrix 
containing periodically distibuted poroelastic inclusions, with low permeabil- 
ities and with imperfect contact on the interface. The micro-models are based 
on Biot's system for consolidation processes in each phase, with interfacial 
barrier formulation. Using the two-scale convergence technique, it is shown 
that the derived system is a general model of that proposed by Aifantis, plus 
an extra memory term. 



1. Introduction 

The interaction between fluid flow and solid deformation in porous media is of 
great importance in petroleum engineering and geomechanics, biosciences, chemical 
processes and many industrial applications [12l HU ED] . 

Some type of porous rocks, like aquifers and petroleum reservoir systems, may 
contain fractures. It is known that flows in such media occur mainly in the fracture 
region and the dominant fluid storage is in the matrix blocks. In this situation, 
reservoir possesses two porous structures, one related to the matrix, and the other 
related to fractures. This notion of double porosity /permeability has first been 
introduced by Barenblatt, Zheltov and Kochina [7] to model the flow of a slightly 
compressible fluid within naturally fractured porous media. The proposed model 
is a system of two partial differential equations in a two-medium description, with 
Darcy's law in each phase, plus exchange extra-terms representing the interfacial 
coupling that results from the interaction, at the micro-scale, between the two 
phases, see (|I.6j) - (|I.7j) below. This was derived under the main assumption that 
the fluid pressure is uniformly distributed at the surface of each phase. 

Generally, fractured rock formations present at the micro-scale high degrees of 
heterogeneity and permeability is mainly determined by the size of the pores and 
the connectedness of the fractures system. So any mathematical modeling of fluid 
flow in such porous media must take into account the rapid spatial variation of 
the phenomenological parameters. Furthermore, from the numerical point of view, 
modeling of such systems at the local scale yields a huge number of discretized 
equations, so computations will be fastidious and intractable. To deal with such 
highly heterogeneous domains, the idea is to replace the medium by an effective 
one. Homogenization techniques, like two-scale convergence method, have been 
used to rigorously derive an effective double-porosity model for the Barenblatt, 
Zheltov and Kochina (BZK) system, see for instance H. Ene and D. Polisevski 
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[14] . However, this model does not take into account the elastic behavior of the 
solid. In fact, a rise in pore pressure of the fluid produces a dilation of the solid 
mass. On the other hand, compression of the medium will increase pore pressure. 
This coupled pressure-deformation was first introduced by Terzaghi [19] in the one- 
dimensional setting and gave the first soil consolidation problem for a homogeneous 
elastic porous medium. Later, M. A. Biot 8. has developed in the multidimensional 
setting a linear theoretical analysis for the behavior of a fluid saturated poroelastic 
medium. The model was based on macroscopic description of the phenomenological 
and physical quantities where the representative volume element is described as 
the superposition of a particle of fluid and a particle of solid. Assuming that 
microstructures are periodically distributed and that the pore scale is very small 
compared to the macroscopic scale, a two-scale asymptotic expansion technique can 
be used to rigorously justify this Biot's model. The microscopic models are based 
on the linear elasticity equations in the skeleton and on the Stokes equations in 
the fluid with appropriate transmission conditions. For more details, we refer the 
reader to the earlier work by Auriault and Sanchez-Palencia [6]. 

Because of the coupling between the deformation and fluid pressure in double 
porosity rocks, which must be understood in order to predict reservoir or aquifer 
behavior, the concept of double porosity has been developed by E.C. Aifantis [T] 
to model oil flow in porous elastic rocks. More precisely, E. C. Aifantis gave a 
phenomenological model for flow of a weakly compressible fluid in a complex and 
heterogeneous medium where a system of partial differential equations is given and 
generalizing Biot's consolidation model by taking into account the basic physics of 
flow through fractured media with interscale couplings. The proposed model reads 
as follows: 

- /iAu - (A + /i)V(divu) + aiVpi + a 2 Vp 2 = f, (1.1) 
cidtpi + aidiv(<9 t u) - i^Api + g(pi - p 2 ) = hi, (1.2) 
cid t p 2 + a 2 div (d t u) - K 2 Ap 2 - g(p x - p 2 ) = h 2 (1.3) 

where u is the displacement of the medium; A and [i are the dilation and shear 
moduli of elasticity, respectively; pi is the pressure of the fluid in phase (z); Cj 
the compressibility, Ki the permeability and an is the Biot- Willis parameters [9]. 
We note that if we let the volume of fissures shrinks to zero so that c 2 ,a 2 ,K 2 ,g 
become negligible then the system (|l.ip - (|1.3|) reduces to the classical Biot system 
with single porosity [8]: 

- /iAu- (A + /i)V(divu) + aiVpi = f, (1.4) 
cidtpi + aidiv (<9 t u) - K x Api = hi. (1.5) 

On the other hand, by neglecting the deformation effects A, /z and ai the system 
(HH])-([I3]) reduces to the BZK model [7J: 

cidtpi - KiApi + g (pi - p 2 ) = hi, (1.6) 
c 2 d t p 2 - K 2 Ap 2 ~ g (pi - p 2 ) = h 2 (1.7) 

Aifantis' theory of consolidation with the concept of double porosity unify then 
the proposed models (|1.4[) - (|1.5[) of Biot for consolidation of deformable porous me- 
dia with single porosity and ()1.6j) - (jl.7p of BZK model for fluid flow through un- 
deformable porous media with double porosity. Note also that a mathematical 
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justification of the Aifantis model has been established in [2]. More precisely, it 
is considered micro-models with periodically distributed poroelastic inclusions, em- 
bedded in an extra poroelastic matrix, with imperfect contact on the interface. The 
micro-model is based on Biot's system for consolidation processes with interfacial 
barrier formulation. The macro-model is then derived by means of the two-scale 
convergence technique and it reads as follows: 



where a, a.i,/3i and Ki are some effective tensors, i = 1,2. See [2] for more details. 
It is then worth pointing out that the Aifantis model (jl.l|l - (jl.3|l can be seen as 
a special case of the homogenized model (|1.8I) - (|1.10I) (ft = a>i = 7J3, 7, being a 
scalar and I3 the identity matrix). 

In this paper, we consider a family of microscopic models for the fluid flow in a 
periodic poroelastic medium made of two constituents : the matrix and the inclu- 
sions, where the material properties change rapidly on a small scale characterized 
by a parameter e representing the periodicity of the medium. We shall make the 
essential assumption that these inclusions have sizes large enough compared with 
the sizes of pores so that it makes sense to consider these media as poroelastic 
materials. 

An interesting question is to investigate the limiting behavior of such media 
when the flow in the inclusions presents very high frequency spatial variations 
as a result of a relatively very low permeability when comparing to the matrix 
permeability, since pore flow velocities in the porous matrix can be high compared 
to movement through the interconnected pore spaces in the inclusions. This leads 
especially to rescale the flow potential in the inclusions by e 2 , as in Arbogast & 
al. [5]. The main objective of this paper is to derive a general model from the 
point of view of homogenization theory. It will be seen that the macro-model is in 
some sense the limit of a family of periodic micro-models in which the size of the 
periodicity approach zero. It is shown that the overall behavior of fluid flow in such 
heterogeneous media with low permeability at the micro-scale may present memory 
terms. It is also shown that in such anisotropic media, with different coupling 
interaction properties in different directions, the Biot- Willis parameters are, as in 
[2], matrices and no longer scalars, as it is usually considered in the poroelasticity 
litterature, since it is assumed there that the medium is homogeneous and isotropic. 

The paper is organized as follows. In the next section [21 we give the geometrical 
setting, the family of the periodic micro-models, and state the main result of the 
paper. Section [3] is devoted to the proof of the main result with the help of the 
two-scale convergence technique. Finally, in section 21 we conclude the paper with 
some remarks. 



The aim of this section is to provide a detailed set up of the studied microstruc- 
ture problem, introduce some necessary notations, basic mathematical tools as well 
as the notion of two-scale convergence, auxiliary problems, and then formulate the 
main result of the paper. 



—diver (u) + aiVpi + a 2 \7p 2 = f , 
d t {cipx +/3i : e (u)) - div (KiVpi) + g(pi - p 2 ) = hi 
d t (c 2 p 2 +Pi-e (u)) - div (K 2 Vp 2 ) - g {pi - p 2 ) = h 2 



(1.8) 
(1.9) 
(1.10) 



2. Setting of the micro-model and main result 
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We consider a bounded and smooth domain in R 3 , e > a sufficiently small 
parameter (e < 1) and Y =]0, 1[ 3 the generic cell of periodicity. We assume that Y 
is divided as Y = Y\ U Y% U T where Y\ , Yi are two connected open subsets of Y and 
r is a smooth surface that separates them. They are such that Y% C Y, Y\ Dia = $j 
r = FT n T2 = (9^2 and dFi = T U 5y. We denote n = ( 

^i)i<i<3 the unit normal 

vector on Y pointing outward with respect to Yf. Let xi>X2 denote respectively 
the characteristic function of Yi, Y% extended by F-periodicity to M 3 . Denote for 

x E 0, x\ (?) = Xi( x / £ ) and set 

Of = {x E O : Xi(x) = 1} and r £ = Of n Of. 

Let = U ee z 3 (^i + e ); * = 1) 2. As in [3], we shall assume that the subset Z\ is 
smooth and connected open subset of M 3 . 

With the above assumptions, the material occupying the domain 0| is then 
embedded in the material occupying Of , and the interface T £ is the boundary dO|. 
We observe that the boundary of Of consists of two parts the outer boundary 
<90 and L £ . Usually, the material Of is referred to the matrix material while the 
material 0| to the inclusions. Note that no connectedness assumption is made on 
the material part 0|. 

Let T > and t E [0, T] denote the time variable. We set the space-time domains 
Q = (0, T) x O, S = (0, T) x T, Q\ = (0, T) x fif , and E £ = (0, T) x T e . 

Let us assume that each phase (Of, 0|) is occupied by a porous and deformable 
material through which a slightly compressible and viscous fluid flow diffuses. Let 
u? denote the displacement of the medium Of, i = 1,2. The equation of motion in 
Of U 0| is given by 

- divaf = f x , in Of, (2.1) 
-divcrj = f 2 , in Of (2.2) 

where of is the stress tensor which satisfies a constitutive equation of linear poroe- 
lasticity of the form |12j : 

of = Afe(uD-afpfI 3 , m Of (2.3) 

and fi E L 2 (O) 3 is the volume distributed force in the corresponding medium, 
i = 1,2. It is assumed that is independent of e. In (|2.3[) . Af and Af are fourth 
rank elasticity tensors, e (•) is the linearized strain tensor, I3 is the identity matrix, 
p e is the pressure and a\ is the Biot- Willis parameter in the poroelastic material 

nf El- 

Let cf (resp. cf) and Kf (resp. .fTf) denote respectively the porosity and the 
permeability of the medium Of (resp. Of). The equation for mass conservation in 
each phase reads as follows: 

d t (c\p\ + a\ divuf ) - div {K{Vp\) = 0, in Of, (2.4) 
d t (cfpf + afdivuf ) - div (iff Vpf ) = in Of. (2.5) 

On the interface T e , we associate to ()2.1j) - (|2.2[) the following transmission condi- 
tions: 

uf = u|, af-n e = af-n E (2.6) 
and to (|2.4[) - (|2.5[) the well-known open-pore conditions: 



(KtVp\) ■ n £ = (iff Vpf ) ■ n £ , (K{Vp\) ■ n £ = -«f (pf - pf ) . (2.7) 
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where n e stands for the unit normal vector on T e pointing outward with respect to 
Of , and g £ is the hydraulic permeability of the thin layer T e . The interface condi- 
tions (I2.7|) are also known as Deresiewicz-Skalak boundary conditions [T3]- Taking 
the limit on the thickness of the thin layer, one can prove that these conditions 
are the only ones that are fully consistent with the validity of the poroelasticity's 
equations throughout heterogeneous media presenting interfaces across which the 
diffusion is discontinuous, see |15) . Observe that when g 6 = oo, (|2.7[) reduces to 
the standard transmission condition, that is a perfect hydraulic contact on the in- 
terface, and when g e = 0, condition (|2.7I) implies no motion of the fluid relative 
to the solid. Here, in this paper we shall assume that none of these conditions are 
fulfilled. See assumption (H4) below. 

On the exterior boundary dil\T £ , we assume the homogeneous Dirichlet bound- 
ary condition: 

ul = and p\ = 0. (2.8) 

Finally, the system (|2.4j) - (|2.8j) is supplemented by the following initial conditions: 

uf(0,.) = 0, p?(0,-)=0mfi?, (2.9) 

ul(0,-) = 0, p§ (0,0=0 hi Jl§. (2.10) 

To deal with periodic homogenization with microstructures, we shall assume the 
followings: 

(HI) There exists y-periodic, fourth rank tensor-valued functions A^ (y), i = 1, 2 
and continuous on R 3 such that 

Af (x) = At ( , x G o, 

and 

(A l (y)E:E)>C(E:E). 
for all y G R 3 and E G .M 3 * 3 (R); 
(H2) There exist ^-periodic real-valued functions Ci (y), i — 1,2 and continuous 
on R 3 such that 

c\ (x) = Ci (—^ , x G fl 

and c t (y) > C > for all y G R 3 ; 
(H3) There exist l"-periodic matrix-valued functions Ki (y), i = 1, 2, continuous 
on R 3 such that 

Kf(x) = K 1 (^j, Kl(x)=s 2 K 2 (^), xen (2.11) 

and 

(K&t) >C|^| 2 , i=l,2 
for all y 6 M 3 and (el 3 ; 
(H4) There exists a function g G C (R 3 ) , ^-periodic such that 

g e (ac) = eg (x/e) , iel 3 and inf g (y) > C > 0. 

t/GR 3 

(H5) The Biot- Willis parameter af is defined a.e. in f2 as follows: 

a\ (x) — a\ for x G Of and (x) = e«2 for x G 0| (2-12) 
where ctj is a positive constant, i = 1,2. 
Here and thoughout this paper, the quantity C denotes various positive constants 
independent of e > 0, of the subscript i = 1, 2 and the microscopic variable y G R 3 . 
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Remark 2.1. We have chosen a particular scaling of the permeability coefficients 
in h2.11\) . This means that the permeability is much larger in the network of 
pores of the inclusions than in the porous matrix rocks. This gives that both terms 
In' \^Pi 1 2 and e 2 J^ e |Vp|| 2 dx have the same order of magnitude and thus lead- 
ing to a balance in potential energies. For more details, we refer the reader to 
Arbogast, Douglas, and Hornung [5] (see also Allaire [3],). In the same way, we also 
have taken a special scaling factor of the Biot- Willis parameters in i 2. 12\) leading to 
a balance in compressibility /dilation between the matrix and inclusions. 

To set the mathematical framework of our Porblem, we need to introduce the 
following spaces: 



2) > 



h = Hi (ny , u = i? (Of) x i? (n: 

£f = {qeH 1 ; q {r = 0} , £| = H 1 (Q s 2 ) , £ £ = £f x £%. 

The space H is equipped with the standard norm: ||v||h = ||e (v) | Iz^ 2 (f2) 3x3 an< ^ & e 
with 

= IIVgiH^n;) +£ 2 HVgallia^) +e ll?i ~ 8a|lz»(r-) ■ 

See S. Monsurro [16]. For a.e. (t,x) G Q, we denote 

u e (t,x) = x\{x)u\{t,x)+x e 2{x)n e 2 {t,x), 
K £ (x) = X \{x)k\{x)+xl{x)M 2 {x) 1 
F(x) = xl(x)h(x) + xl(x)h(x). 

Note that, thanks to the transmission condition (|2.6p , the displacement u £ (t, •) lies 

in H for a.e. t e (0,T). 

Throughout the paper the following notation will be used: if F is any functional 

space then I? T (F) denotes the Bochner vector-valued functions space defined by 

If T (F) = LP (0, T; F) 

The weak formulation of (I2.4[) - (|2.10|) can be read as follows: 

Find (u £ ,p e ) G Lf (H) x L\ (£ e ), such that p £ = (pf G Lf (L £ ) , 

d t {c\p\ + axdivu 5 ) G L% {El*) , d t (c e 2 p E 2 + ea 2 divu e ) G L% (£ £ 2 *) 

and for all v G H, (91,(72) £ £ £ , we have 

/ A £ e (u £ ) e (v) dx + / aiVpfvdx + / a £ 2 Vp £ 2 vdx = / i £ vdx, (2.13) 
Jn Jn\ Jn* Jn 

(d t (c\p\ + aidivu e ) , qi) S s. £c + / K £ Vp\Vqidx + 

1 ' 1 Jn\ 

(d t (c 2 p 2 + £a 2 divu e ) , 52)5*. £ e + K 2 \7p 2 \7q 2 dx + 

2 ' 2 Jn* 

I g £ (p\-p e 2 ){qi-q2)ds £ {x)=Q, (2.14) 

JV* 

u £ (0, •) = 0, xi (-)P\ (0, •) + X2 (-)P £ 2 (0, •) = a.e. in fl (2.15) 

Here and throughout this paper dx and ds e (x) stand respectively for the Lebesgue 
measure on R 3 and the Hausdorff measure on T e . 
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Using assumptions (H1)-(H5), we can establish the following existence and unique- 
ness result whose proof is a slight modification of that given by R. E. Showalter 
and B. Momken [TS] and therefore will be omitted. 

Theorem 2.2. Assume that (H1)-(H5) hold. Then, for any sufficiently small e > 
and f £ £ L 2 (Q), there exists a unique couple (u e ,p e ) G L^? (H) x L\ (£ e ), solution 
of the weak system 12.13\) - ![2.15\) . such that 

l!u £ IU»( H ) + Wfhusn + Iblw) < o. (2.16) 

Now, thanks to the a priori estimates ([2.16[) , one is led to study the limiting 
behavior of the sequence (u £ ,p £ ) as e approaches 0. To do this, we shall use the 
two-scale convergence technique that we shall recall hereafter. 

First, we define C#(Y) to be the space of all continuous functions on M 3 which 
are F-periodic. Let the space L 2 ^ (Y) (resp. L 2 ^ (Yi), i = 1,2) to be all functions 
belonging to Lf oc (M. 3 ) (resp. L 2 oc (Zi)) which are F-periodic, and (Y) (resp. 
Hg (Yi)) to be the space of those functions together with their derivatives belonging 

toZ|(Y) ( res P- L % ( Z i))- 

Now, we recall the definition and main results concerning the method of two-scale 
convergence. For more details, we refer the reader to [51 141 117], 

Definition 2.3. A sequence (v £ ) in L 2 (f2) two-scale converges to v G L 2 (17 x Y) 

(we write v e v) if, for any admissible test function tp G L 2 (Q;C#(Y)), 



lim 



i / v £ (x) tp (x, - J dx = / v (x,y) ip (x,y) dxdy. 
°Jn \ eJ J QxY 

Theorem 2.4. Let (v £ ) be a sequence of functions in L 2 (Q) which is uniformly 
bounded. Then, there exist v G L 2 (Jl x Y) and a subsequence of (v £ ) which two- 
scale converges to v. 

Theorem 2.5. Let (v £ ) be a uniformly bounded sequence in iJ 1 (J7) (resp. Hq(Q)). 
Then there exist v G H 1 (17) (resp. H^(Q)) and v G L 2 (17; H^(Y)/R) such that, up 
to a subsequence, 

v £ 2 ^ v; Vv £ 2 ^ s Vv + V y v. 

Here and in the sequel the subscript y on a differential operator denotes the 
derivative with respect to y. 

Theorem 2.6. Let (v £ ) be a sequence of functions in H 1 ^!) such that 

\\v e \\mn)+z\\Vv £ \\ L2{n)3 <C. 

Then, there exist v G L 2 \Q; Hg(Y)\ and a subsequence of (v £ ), still denoted by 
(v £ ) such that 

v £ v, eVv £ 2 — ^ V y v 
and for every tp G £> (17; C#(Y)), we have: 



lim 



i/ v £ (x)<p (x, -) ds £ (x) = / v(x,y)(p(x,y)dxds(y) 



Here and in the sequel ds (y) denotes the Hausdorff measure on T. 
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The notion of two-scale convergence can easily be generalized to time-dependent 
functions without affecting the results stated above. According to [10], we give the 
following: 

Definition 2.7. We say that a sequence (v £ ) in L 2 (Q) two-scale converges to 
v G L 2 (Q x Y) (we always write v e 2 -^ s v) if, for any if G L 2 {Q;C#{Y)): 



(t,x) ip (t,x, — J dtdx = I v (t, x,y) ip (t, x,y) dtdxdy. 



lim 

' *°JQ v c ' JQxY 

Remark 2.8. The results stated above still hold for the case of time- dependent 
sequences. For if (v £ ) is a uniformly bounded sequence in L 2 {Q), there exists then 

v G L 2 (Q) such that, up to a subsequence, v e 2 — ^ v in the sense of Def. \2. 7[ 
Moreover, if (v £ ) is uniformly bounded in L|, (0)), then up to a subsequence, 

there exist v G L\ {H 1 (fi)) and v a G L 2 (Q; H# (Y) /Rj such that v 



2-s 

— 1 I! ano 



Vu e 2 — ^ Vu + VyVQ. On the other hand, if a sequence (v s ) is such that 

IKIIl2 ( q)+eIIW|| l2(q) <c, 

then, up to a subsequence, there exists v G L\ 

(## { Y )j such that 

v E v and 

eV y v e 2 — * VyV. 

In order to state the main result, we shall give in the sequel some notations. 
Let us first introduce the three following auxilliary problems. For j,k G {1,2,3}, 

let w jfc G (^H# (Y) /M^ be the solution to the following microscopic system: 

-div a (Axe y {w^ k + d jk )) = a.e. in Y u 
-div y (A 2 e y (w^ fe + d jk )) = a.e. in Y 2 , 
Aie y (w J& + d jk ) ■ n = A 2 e y (w jk + d? k ) ■ n a.e. on T, 
y i — > w J Y"-periodic 

where d jk (y) = {yjO~ik) l<l<5 and {5k j) is the Kronecker symbol. For j = 1, 2, 3, let 
7Tj G H 1 {Yi) /R be the solution of the following stationary micro-pressure equation: 

-dWy {K x (VtTj + ej)) = in Yi, 
K\ (VtTj + ej) ■ n = on T, 
y i — > Tij Y-periodic 

where e.j is the j th vector of the canonical basis of R 3 . Let £ G L^? {YiYj be 
the unique solution to the following non micro-pressure problem of the Robin type: 

{ d t (c 2 C) - divy {K 2 VyQ = a.e.in (0, T) x Y 2 , 
K 2 V y C • n= - g (y) [1 - C] a.e. on E, 
y i — > £ Y"-periodic, 
[ C(0,y) = 0, a.e. y G Y 2 . 

Now, let us define the homogenized fourth rank tensor A = fa,-, ,-,,-,,-.),^ .... 
where the coefficients are given by 

3 f 

ajihhu = / a hhkik 2 (y) (5 jl k 1 5 32 k 2 + e kl k 2 , y (w J « 4 ) (y)) dy. 

k 1 ,k 2 =i 
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Here (ajkim) arc the coefficients of the elasticity tensor A which are given by 

A (y) = xi (y) Ai (») + X 2 (y) A 2 (y) (2.17) 

for a.e. y G Y, and e^y (•) is the linearized elasticity strain tensor where the 
derivatives are taken with respect to the microscopic variable y. Let also define the 
following homogenized tensors: 

er(u) = (^ fc (u)) , K = (R jk ) , B = (b jk ) , A = (X jk ) (2.18) 

where for j,k S {1, 2, 3} 

3 

5j'fc(u) = a J fc im e im (u), (2.19) 

Z,m=l 



= / Ki (y) (V„7r,- + e 3 ) (V7T fe + e k ) dy, (2.20) 
b jk = ai f + y ^fc (») njds (y)^ , (2.21) 



' 1 i / > ; | 6 jl 6 k i + -Q^-\dy. (2.22) 



Here I Y I denotes the volume of Y and ( ) are the components of w u . 

V /i<;<3 

Finally let us define the following averaging quantities 

f = |Yi|fi + |Y 2 |f 2 , (2.23) 
ci(y)dy, (2.24) 



IY 1 

9 = j g{y)ds{ y ) (2.25) 
and the time-dependent functions 

0(t,r) = a 2 J d t ((t-T,y)nds(y), (2.26) 

r?(t,r) = - [ g(y)d t ((t-T,y)d S (y). (2.27) 



With these notations, we are now ready to give the main result of the paper: 

Theorem 2.9. Let(u £ ,p £ ) 6 i°° (0, T; H) xi 2 (0, T; £ e ) be the solution of the weak 
system (2.13$ . Then, up to a subsequence, there exists (u,f>) € I? (0, T; Hj (f2) x i?Q (f2)) 

u £ urn L 2 (0, T; iJp 1 (ft)) weaifcZj/, 

p\ — 1 pi m i 2 (Q) weakly, 

p\ P2 (y) dy in L 2 (Q) weakly, 

Jy 2 

where p = (pijy 2 P2 (y) dy^j , 

p 2 (t,x,y) = / pi(T,x)d t C(t-T,y)dT, a.e. (t,x,y) £ Q x Y 2 . 
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and the couple (u,pi) is a solution to the homogenized model: 

— diva(u) + BVpi + / 8 (t, t) pi (r, x) dr= f , a. e. in Q, 
Jo 

dt (cpi + A : e (u)) — div (^KVpij + gp\ — J n (t,r) p\{r,x) dr=0, a.e. inQ, 

u= and ICVpi ■ v=0 a.e. on E, 

u (0, x) = a.e. in £1, p\ (0, x) = a.e. in Q, 

Here a, B, 9, f, c, A, K, g and n are given in 12.18]) - f2.27\ ). 

3. Proof of the main result 

As a direct application of the theorems listed above (Thms 12.4112.61) and the a 
priori estimates (|2.16[) . we give without proof the following two-scale convergence 
result concerning the solutions (u £ ,p £ ) of the Problem (|2.13p - (|2.15|) . 

Theorem 3.1. There exists a subsequence of (u £ ,p £ ), solution of k2.13\) - l2.15)) . 
still denoted (u £ ,p £ ), and there exist 

ue^(H), ue^(L 2 (J];4(r)/i)) 3 
pi g l§? {Hi (n)) , n g £ 2 (Q; (r) /k) 

and 

Kel? (L 2 
swc/i i/ia£, /or a.e. i G (0,T), 

u e (t,-) 2 -^u (*,•), (3.1) 

x!p! (t, •) 2 -^ s XlPi (*, •) . (3-2) 
X?Pl (*, •) ^ X2P2 (t, •) (3.3) 



in the sense of Def. \2.3[ and 

<9u £ 2-s du du „ „ „ ,„ ,. 

cxj cxj ay., 

xfVp? 2 - s xi (Vpi + Vj.pi) , (3.5) 

exl^pl X 2^vP2 (3.6) 
in the sense of Def. \2. 7[ Moreover, the following convergence holds: 

lim / s {p{ - p £ 2 ) ip £ dtds 6 = / (p 1 -p 2 )ipdtdxds, (3.7) 

for any ^eD (Q; C# (Y)) with %p £ (t, x) = tp (t, x, x/e). 

To determine the limiting equations of the system (|2.13l) - (|2.15l) . we begin by 
choosing the adequate admissible test functions. Let v 6 (x) = v (x) + ev \ x, —^j 

where v G V (ft) 3 and v G V (Q;e^ (Y)^ 3 . Let also qf (t, x) = ip x (t,x)+e<p 1 (t,x, 
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andg|(t,x) = ip 2 (t,x,-J where cpi e 2?((0,T) x Cl) and ip 2 , <p\£V (Q;Cg (y)). 
Taking v = v £ in (|2.13|l . we have: 

/ f £ v £ dx = / A £ (x)e(u £ )e(v E )dx + / a{sjp\v e dx + e / a 2 Vp e 2 v e dx 
Jn Jn Ja\ Jn* 

= J A £ (x) e (u £ ) (e (v) (x) + e y (v) (x, dx + 

[ {a l xl{x)V P l+ea 2 x e 2 {x)V P e 2 )^{x)dx + eR\, (3.8) 
Jn 



where 



R\ = j A £ (x) c (u £ ) c x (w) (x, — ^ + ai ^ x! ( x ) Vpfw ^r, — ^ da; 
+£"2 / X2 ( x ) Vp^w (x, da. 



Observe that i?f = O (1). 

Now, we pass to the limit in (13.8[) . In view of (13.41) . and since A* (e (v) + e y (v)) 
is an admissible test function, the first integral in the l.h.s. of (|3.8[> converges to 



A (c (u) + c y (u)) (e (v) + e y (v)) dxdy (3.9) 

where the tensor A (y) is given by (|2.17|) . In view of Divergence Lemma and (|3.5jl - 
(|3.6p . the second integral of the l.h.s. of Q3.8P tends to 

ai j (Vpi + VyPi) v {x) dxdy + a% / V y p 2 v (x) dxdy 

= ai\Yi\ / Vpiv (x) dx + / (aipi + a 2 p 2 ) (v ■ n) cfeds, (3.10) 

By Theorem 12.41 it follows that 

lim / f £ v £ (s) = lira f / f £ Cx) w (x) dx + e [ f £ (s) v far, da; | 
£ ^° Jn £ ^° V^n v £ y / 

= / fv(ar)da; (3.11) 

where f is given by (|2.23[) . Thus, collecting these limits (|3.9[) - (|3.1ip . we obtain the 
limiting equation of 



/ A [e (u) + c y (u)] [e (v) + c y (v)] dxdy + an \Yi\ / Vpivdx 
JnxY Jn 

+ / (aipi + a 2 p 2 ) (v • n) dxds = I fvdx (3-12) 
Jnxr Jn 

which is valid for a.e. t £ (0, T). Next, we proceed to get the limiting equation of 
(12.141) . Taking q\ = 1\ an d q 2 = q 2 in (|2.14[) . integrating by parts over (0,T) and 
taking into account the initial conditions (|2.15[) . we obtain 



(c? {x)p\ + aidivu £ ) d t <fii (t, x) dtdx — / cf> (x)p 2 d t ip 2 (t, x, — J dtdx + 

JQI \ eJ 

J K\ ( Vp £ (v^i (t, s) + V y ip! (t, x, ) dtdx + 



12 ABDELHAMID AINOUZ 



e 

where 



ek 2 \ Vp|V tf V72 (t, x, dtdx + 
5 (J) ifi ~ Pi) (vi (*. x) - <P2 (t, x, ~) ) dids E + ei?2 (3. 



— {c\ {x) p\ + aidivu 6 ) dt^i (t, x, — dtdx + 
/ — a2divu e d t <p2 (t, x, — ^ dtdx + 
J K Y (^J Vp\V x (pi (t, x, ^ dtdx + 
e J K 2 (^J Vp e 2 V x tp 2 (t,x, dtdx + 
9 (^) (Pi -PD'PI (t,x)< 



i dtds e 



The first integral of Q3.13P is equal to 



J ~ Xl (f) ( Cl (f ) p i + aidivu e ^ d t (pi (t, x) dtdx, 
and thanks to (|3.2p and (|3.4p . converges to 

"Xi (y) (ci (y)pi +ai (divu + div^u))^^! (t 7 x) dtdxdy. 



'QxY 

In a similar way, by (|3.3[) and (|3.4p . it follows that 



(x) p £ 2 dtf2 (t, x, - ) dtdx -> / X2 (?/) c 2 (y) £>2<9t<P2 (*, y) dtdxdy 



e 

Now, in view of p. 51) one can deduce that 



dtdx 



y Xl ( #1 (J) Vpf (V<y9l (t, S) + Vy(f>l (t, X, J) ) dtdx -> 

/ Xi (v) K i (v) (Vpi + Vj,pi) (V<^ (f, x) + V y £i (t, x, y)) didxdy 

JQxY 



IQxY 

and thanks to (|3.6p . we also get 



efc 2 Vp|Vj/<^2 (*,x, didx 
*2 #2 eVp e 2 V y <p 2 (t, x, ^ d^dx 
X2 (y) i^2 (y) Vp 2 V a (p 2 (*, as, y) dtdxdy. 



}xY 

By virtue of (|3.7p . we find that 

x 
e 



9 (-) (Pi - Pa) (^1 (*> x ) " ^2 (<, x, ■ 



dtds e -)• 
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9 (y) {pi ~ P2) ((pi (f, a;) - (fi2 (t, x, y)) dtdsdy. 

Qxr 

As before, we observe that R 2 = O (1) and, by collecting all the preceeding 
limits, we get the following limiting equation of (12. 14|) : 

— (ci (y) pi + ai (divu + div^u)) dt<pidtdxdy + 
/ Ki (y) (Vpi + Vypi) (V<pi + Vyfii) dtdxdy + 

JQxY! 

(-c 2 (y) p 2 d t (p 2 + K 2 (y) V y p 2 V y <p 2 ) dtdxdy + 

9 (y) (Pi - P2) (<fi - P2) dtdsdy = 0. (3.14) 

By density argument, the equations (|3.12l) and (13. 14|) still hold true for any 
(v, v) e Hxl 2 (ft, H 1 (Y) /R) 3 and (<p u y x , <p 2 ) G L% (H 1 (Q)) xL 2 (Q; JT* (Y) /m) x 

i 2 (<3; (y)) • We can summarize the preceding by observing that these equa- 
tions are a weak formulation associated to the two-scale homogenized system (|3.15[1 - 
(|3.31|) . Indeed, integrating by parts in (|3.12[) and p.!4[) . we obtain the following 
system: 

~dW y (Ai [e (u) + e v (u)]) = a.e.in QxY 1: (3.15) 
-div y (A 2 [e (u) + e y (u)]) = a.e.in Q x Y 2 , (3.16) 

-div (^J A [e (u) + c y (u)] dyj + a>i | Yi | Vpi + 

y (aipi + a 2 p 2 ) nds= f a.e. in Q, (3-17) 

and 

-div a (Kx (Vpi + V y pi)) = a.e. in Q x Yi, (3.18) 

3 t (c 2 p 2 ) - div y (K 2 \7 yP2 ) = a.e.in Q x Y 2 , (3.19) 

J 9 (y) [Pi - Pi] (y) = a.e. in Q, (3.20) 

with the transmission and boundary conditions: 

Ai [e (u) + e y (u)] ■ n= A 2 [e (u) + e y (u)] ■ n a.e. on Q x T, (3.21) 



(Ki (Vpi + V y pi)) • n= a.e. on Q x T, 



(3.22) 
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(Ki (Vpi + V y pi)) ■ v= a.e. on (0, T) X dtt X Fi, (3.23) 

K 2 V y p 2 ■ n= - g (y) [pi - p 2 ] a.e. on Q x T, (3.24) 

u= a.e. on 90, (3.25) 

y — ► u, pi,P2 F-periodic, (3.26) 
and the initial conditions: 

u(0,x)= a.e. in Q, (3.27) 

u(0,x,y) = a.e. in Q x F, (3.28) 

Pi (0,x)= a.e. in 0, (3.29) 

pi (0, a, y) = a.e. m Q X 7i (3.30) 

p 2 (0,i,y)=0a.c. inOxF 2 . (3.31) 

Now we decouple the system (|3.15[) - (|3.31[) . In view of the linearity of the two first 
equations (|3.15|) - (|3.16|) . we can write that, up to an additive constant: 

3 

u(i, x, y) = *j (u) (t, x) w« (y) + C te , a.e. (t, x, y) G Q x F, (3.32) 

i,j=l 



where, for i,j G {1,2,3}, w 4 - 7 G ( i/^ (F) /KJ is the solution to the following 
microscopic system: 

-div a (Aie„ (w ij + d 4J )) = a.e. in Y 1 , 
-dw y (A 2 e y (-w ij + d l] )) = a.e. in Y 2 , 
Aie y (w y + d ij ) ■ n = A 2 c y (w' y + d ij ) ■ n a.e. on T, 

y i — > w IJ F-periodic. 

Here d fcz = (yK^u) 1<i<3 and is the Kronecker symbol. 

Similarly, in view of (|3.18p . (I3.22[) and (|3 . 26[) one can write that: 

3 q 

p l (t,x,y)=^2-£-{t,x)T H (y) + C te , a.e. (t, x, y) G Q x Y u (3.33) 

2—1 ? 

where, for i = 1,2,3, the micro-pressure 7Tj G -ff 1 (Yi) /R is the solution of the 
following stationary equation: 

-div a (ifi (Vvr 4 + e,)) = in Y u 
Ki (V7Tj + ej) • n = on T, 
y i — > -Ki F-periodic. 

Here is the i th vector of the canonical basis of M 3 . Let us denote 

A = (Si 1 i 2 t3i 4 )i<j li j 2i i 3i i 4 <3 i 
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3 r 

/ a ^23U2 (y) {^h § i2h + e hh,v (w l3U ) (y)) dy, 

A -I JY 



31,32 = 1 ' 

where (aiji m ) are the coefficients of the elasticity tensor A and 
, . 1 / dwi dvjj \ . 

e ^ (w) = H% + ^-J' w=K)i ^ 3 - 

Let also define the effective stress tensor 

3 

o"( u ) = (M u ))i<ij< 3 ' °ij( n ) = Y a«ime; m (u), 

/,m— 1 

the effective permeability tensor 



K = 



(k l3 ) , K i:j = I Ki (y) (Vym + ei) (Vttj + e s ) dy, 



l<ij<3 J Yl 

the effective Biot- Willis matrices: 

B = (bij) , hj = ai Mli | Sij + TTj 



3 (y)nids(y) ) , n = (m)^^ 

A = (K 0l<i,j<3 ' X *3 = "1 / H ( ^"A™ + ) 



and finally the following averaging quantities 



ci (y) dj/, 9= 9 (y) ds (y) . 
Then from (|3.32t - (|3.33p we deduce the homogenized system 



-div<r(u) + BVpi + a 2 J p 2 nds (y) — f a.e. in Q, (3.34) 
dt (cpi + A : e (u)) - div (^KVp-^j + gpi 

g (y) p 2 ds (y) = 0, a.e. in Q, (3.35) 



d t (c 2 p 2 ) - div y (K 2 V yP2 ) = a.e. in Q x Y 2 , (3.36) 
c 2 V v p 2 n= - g (y) [p x - p 2 ] a.e. on Q x T, (3.37) 
u= and KV Pl ■ v = a.e. on (0,T) x S, (3.38) 
y i — > p 2 F-periodic, (3.39) 
u (0, x) = a.e. in O, pi (0, i) = a.e. in ft, (3.40) 
p 2 (0,x,y)= a.e. in ft x F 2 . (3.41) 
Now, we are going to establish a relation between the two pressures p\ and p 2 - 
To this aim, let ( e i°° (o,T;iJ^(F 2 )) be the unique solution to the following 
microscopic and non homogeneous Robin problem: 

d t (caC) - div y (K 2 V a C) = a.e.in (0,T) x Y 2 , 
K 2 V y C ■ n= - g (y) [1 - C] a.e. on E, 
y i — > £ F-periodic, 
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C(0,y) = 0, a.e. yeY 2 . 

Since C2,K2,g are time-independent and p\ is independent of y, using the Laplace 
transform method, one can then easily see that 

pa(t,x,y)= / pi(T,x)d t ((t-T,y)d,T, a.e. (t, x, y) e Q x Y 2 . (3.42) 
Jo 

Therefore, the homogenized system (|3.34[) - (I3.41[) can be rewritten as 

— divcr(u) + BVpi + I 8 (t, r)pi (t, x) dr— f a.e. in Q, 
Jo 

dt (cpi + A : c (u)) - div LfiTVpiJ + gpi 
t 

T)(t,T)pi(r,x)dT=0, a.e. in Q, 

u= and KVp! ■ v= a.e. on (0, T) x dQ, 
u (0, x) — 0, pi (0, x) = a.e. in ft, 

where we have denoted 

(t, t) = a 2 J dt( (t - t, y) nds (y) , 

v(t,T) = J g(y)d t ((t-T,y)ds(y). 

Finally, let us observe that the overall pressure of the fluid flow in the microstructure 
model which is given by 

P 6 {t, X ) = X\ {x)p\ (t,x) + X2 {X)P B 2 {t, X) 

for a.e. (t,x) € Q, two-scale converges to Xi {v)Pi (t, %) + X2{y)P2 (t, x, y), and 
thanks to (|3.42p . converges then weakly in L 2 (Q) to 

\Y 1 \p 1 (t,x)+ / p 1 (T,x)d t ({t-T,y)dydT. 
Jo Jy 2 

This concludes the proof of Theorem 12.91 

4. Conclusion 

We have used the homogenization theory to derive a macro-model for fluid flow 
in composite poroelastic with microstructures, in which inclusions are fully embed- 
ded and with very low permeabilities. We have shown that the overall behavior 
of fluid flow in such heterogeneous media with low permeability at the micro-scale 
may present memory terms. We also have shown that in such cases, the Biot- Willis 
parameters are, as in matrices and no longer scalars, as it is usually consid- 
ered in the poroelasticity litterature, since it is assumed there that the medium is 
homogeneous and isotropic. Nevertheless, anisotropic media may present different 
coupling interaction properties in different directions at the micro-scale, and which 
lead at the macro-scale to such anisotropic Biot- Willis parameters. Finally, let us 
mention that the result of the paper remains valid if one considers non homogeneous 
initial conditions or with any volume distributed source densities in each phases. 




HOMOGENIZATION OF A DOUBLE POROSITY MODEL IN DEFORMABLE MEDIA 17 



References 

[1] E.C. Aifantis; On the response of fissured rocks, Develop. Mech. 10 (1979), 249-253. 

[2] A. Ainouz; Homogenized double porosity models for poroelastic media with interfacial flow 

barrier, Mathcmatica Bohcmica 136 (2011), 357-365. 
[3] G. Allaire; Homogenization and two-scale convergence, SIAM J. Math. Anal. 23 (1992), 

1482-1519. 

[4] G. Allaire, A. Damlamian and U. Hornung; Two scale convergence on periodic surfaces and 
applications, in : Proceedings of the International Conference on Mathematical Modelling 
of Flow through Porous Media (May 1995). A. Bourgeat et al. eds., 15-25, World Scientific 
Pub., 1996. 

[5] T. Arbogast, J. Douglas and U. Hornung; Derivation of the double porosity model of single 

phase flow via homogenization theory, SIAM J. Math. Anal. 21 (1990), 823-836. 
[6] J.L. Auriault and E. Sanchcz-Palcncia; Comportement macroscopique d'un Milieu Poreux 

Saturee Deformable, J. Mec. 16 (1977), 575-603. 
[7] G. I. Barenblatt, I. P. Zheltov and I. N. Kochina; Basic concepts in the theory of seepage of 

homogeneous liquids in fissured rocks, Prikl. Mat. Mckh. 24 (1960), 852-864. 
[8] M. A. Biot; General theory of three dimensional consolidation, J. Appl. Phys. 12 (1941), 

155-169. 

[9] M. Biot and D.Willis; The elastic coefficients of the theory of consolidation, J. Appl. Mech. 
24 (1957), 594-601. 

[10] G. W. Clark and R. E. Showalter; Two-scale convergence for a flow in a partially fissured 

medium, Elcctr. J. of Diff. Equations. 1999 (1999), No. 02, 1-20. 
[11] S. C. Cowin; Bone poroelasticity : A Survey Article, Journal of Biomechanics 32 (1999), 

217-238. 

[12] O. Coussy; Poromechanics, John Wiley & Sons Ltd, 2004. 

[13] H. Deresiewicz and R. Skalak; On uniqueness in dynamic poroelasicity, Bull. Scismol Soc. 

Amcr. 53 (1963), 783-788. 
[14] H.I. Ene and D. Polisevski; Model of diffusion in partially fissured media, Z. angew. Math. 

Phys. 53 (2002), 1052-1059. 
[15] B. Gurevich and M.Schoenberg; Interface conditions for Biot's equations of poroelasticity, J. 

acoust. Soc. Am. 105 (1999), 2585-2589. 
[16] S. Monsurro; Homogenization of a two-component composite with interfacial thermal barrier, 

Advances in math. Sci. and Appli. 13 (2003), 43-63. 
[17] G. Nguetseng; A general convergence result for a functional related to the theory of homog- 
enization, SIAM J. Math. Anal., 20 (1989), 608-623. 
[18] R. E. Showalter and B. Momken; Single-phase flow in composite poroelastic media, Math. 

Methods Appl. Sci. 25 (2002), 115-139. 
[19] K. von Tcrzaghi; Die berechnung der durchlassigkeitsziffer des tones aus dem verlauf der 

hydrodynamischen spannungserscheinungen, Sitz. Akad. Wissen., Wien Math. Naturwiss. 

Kl. Abt. Ha 132 (1923), 105-124. 
[20] H. F. Wang; Theory of linear poroelasticity with application to geomechanics and hydrogeol- 

ogy, Princeton University Press, Princeton, 2000. 

Abdelhamid AINOUZ 

Laboratoire AMNEDP, Mathematics Department, USTHB, BP 32 ElAlia, Bab-Ezzouar, 
Algiers, Algeria. 

E-mail address: aainouzOusthb . dz 



